Conduction of Ultracold Fermions Through a Mesoscopic Channel 
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In a mesoscopic conductor electric resistance is detected even if the device is defect-free. We 
engineer and study a cold-atom analog of a mesoscopic conductor. It consists of a narrow channel 
connecting two macroscopic reservoirs of fermions that can be switched from ballistic to diffusive. 
We induce a current through the channel and find ohmic conduction, even for a ballistic channel. 
An analysis of in-situ density distributions shows that in the ballistic case the chemical potential 
drop occurs at the entrance and exit of the channel, revealing the presence of contact resistance. 
In contrast, a diffusive channel with disorder displays a chemical potential drop spread over the 
whole channel. Our approach opens the way towards quantum simulation of mesoscopic devices 
with quantum gases. 



The quantum simulation of models from solid-state 
physics using cold atoms has seen tremendous progress 
over the last decade [U|2]. Still, there are only limited 
analogies to the concept of conduction, which is at the 
core of mesoscopic solid-state physics. To close this gap 
it would be highly desirable to connect a probing region 
in a cold atom experiment to external incoherent reser- 
voirs. This would lead to directed transport, the control 
of which is the basis of electronics. In such an intrin- 
sically open configuration, boundary conditions play a 
crucial role, as in the Landauer theory of transport [3]. 
So far, the transport properties in cold atom systems 
have been investigated by observing the response of the 
system to variations of the external potential [4UT3]. or 
by monitoring the coherent evolution of bimodal Bose- 
Einstein condensates [T4UT6] as a response to a bias. Ex- 
tending the concept of quantum simulation to conduc- 
tion requires the engineering of macroscopic reservoirs, 
an atom battery or capacitor connected to the conduc- 
tor [T7HT9] . Such experiments would not only provide 
cold atoms with a faithful modeling of conduction exper- 
iments in solids, but also open the possibility to study 
phenomena that emerge at the contacts of mesoscopic 
systems [20] . 

We report on the observation of atomic conduction be- 
tween two cold atom reservoirs through a mesoscopic, 
multimode channel. Our observations are three fold: (i) 
After preparing the two reservoirs with different atom 
numbers, we measure the atom current as a function of 
time. The system slowly evolves towards equilibrium, re- 
alizing a quasi-steady state at each point in time. This 
allows us to extract the current as a function of the 
atom number difference, (ii) Using in-situ imaging of 
the atomic density in the channel, we establish that the 
chemical potential drop occurs at the contacts between 
the channel and the reservoirs, as predicted by the Lan- 
dauer theory, (iii) Alternatively, we prepare a strongly 
disordered channel having the same macroscopic conduc- 
tion properties, and observe that the chemical potential 
drop is spread over the whole channel, revealing a finite 
linear resistivity. 



Our measurement is made possible by the separation 
of scales in our trap geometry, as illustrated in figure [I] 
The experimental configuration consists of two identical, 
macroscopic cold atom reservoirs, which contain the ma- 
jority of the atoms and feature fast equilibration dynam- 
ics. They are connected by a multimode channel, which 
contains a negligible fraction of the atoms, and supports 
a few quantum states in the z-direction, while it has the 
same extension as the reservoirs in the x-direction, mak- 
ing it quasi-two-dimensional. A high-resolution micro- 
scope is used to detect and influence the atoms at the 
micrometer scale [21] . 
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FIG. 1: Experimental configuration. A macroscopic 6 Li cloud 
is divided into two reservoirs separated by a narrow channel. 
The channel is imprinted using the two lobes in the intensity 
profiles of a nearly TEMni-mode laser beam at the wavelength 
of 532 nm, created with a holographic plate. The distance 
between the two lobes is 18 /xm, and the waist of the beam 
in the ^/-direction is 30 /im. A microscope objective (NA 
= 0.55) is used to observe or manipulate the atoms in the 
channel. When Nl is larger than Nr, an atomic current / 
flows through the channel (represented by the arrow on the 
figure) . 



We prepare quantum degenerate gases containing iV to t 
= 4 • 10 4 6 Li atoms in each of the two lowest hyper- 
fine states at a temperature of 0.36(18) Tp, where Tp — 
700 nK is the Fermi temperature in a combined optical 
and magnetic trap [21] . A laser beam propagating along 
the x-direction is focused on the center of the atomic 
cloud. The beam has a nodal line in the middle of its 



intensity profile and produces a repulsive potential for 
the atoms, which is tightly confining in the z-direction 
[22| [23] . Oscillation frequencies of up to 3.9 kHz along 
the z-direction are achieved, see figure 1. 

Figure [2]A presents a typical absorption picture of a 
cloud in the presence of the channel. We observe two 
clouds clearly separated by a low density region, reveal- 
ing the presence of the channel and confirming that it 
contains a negligible fraction of the total atom number 
(smaller than 0.01). 




FIG. 2: Atomic reservoirs connected by a mesoscopic channel. 
A : Absorption image of the atoms before the imbalance is 
applied. The image is taken after 1 ms expansion in the x-z 
plane. The dark region at the center reveals the presence of 
the mesoscopic channel. The dashed line tracks the position 
of the channel on both panels. B : Density difference between 
the cloud in the unbalanced configuration, as the current sets 
in, and the equilibrated cloud. Excess of atoms is displayed in 
orange, and lack of atoms is displayed in green. The imbalance 
of the reservoirs AN/Ntot is set to 0.2. 

The conduction measurement proceeds as follows. We 
create an asymmetry in the potential by applying a con- 
stant magnetic field gradient of 2.5 mTm" 1 along the 
y axis. This is done during the evaporation process and 
eventually results in an imbalance AN/N tot ~ 0.2, where 
AN is the number difference between the two reservoirs. 
After evaporation, the confining potential of the trap is 
increased and a uniform magnetic field is set to 47.5 mT. 
At this value, the scattering length of atoms in the two 
internal states is -100 ao, with ao being the Bohr radius. 
This ensures that the collision rate is sufficient to main- 
tain thermal equilibrium in each reservoir on a time scale 
of c± 30 ms. It also ensures that the mean free path 
(~ 1.3 mm) is much larger than the length of the channel, 
making it ballistic. The symmetry of the trapping poten- 
tials is then restored by switching off the magnetic field 
gradient in 50 ms, a time longer than the internal ther- 
malization time of each reservoir, but short compared to 
the timescale of equilibration of the populations of the 
two reservoirs. Figure [2j3 shows the difference between 
an absorption picture taken with and without imposing 
an imbalance. The left reservoir is seen to contain an 
excess of particles compared to the balanced reservoirs 
situation, and the right reservoir shows a deficit of par- 
ticles. 

The equilibrium of the whole system is characterized 
by a balanced population of both reservoirs, thus after 
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FIG. 3: Observation of ohmic conduction. A: measured im- 
balance between the two reservoirs as a function of time. The 
solid line is an exponential fit to the data. B : Current as 
a function of number difference between the two reservoirs, 
measured from the exponential fit of figure A, for two differ- 
ent confinements in the channel. A small offset obtained from 
the fits in figure A, which is due to a slight misalignment of 
the channel with respect to the center of the trap, has been 
substracted. Circles : maximum center frequency along z set 
to 3.9 kHz, triangles : 3.2 kHz. The lines are linear fits to the 
data. 



restoring the symmetry of the trap an atomic current sets 
in through the channel. Figure [3]A^ presents the time evo- 
lution of AN/Ntot, with the oscillation frequency along 
z in the channel set to 3.9 kHz. We observe an expo- 
nential decay (solid line on figure [3^) , with a time con- 
stant of 170(14) ms. This exponential shape suggests a 
direct analogy with the discharge of a capacitor through 
a resistance. Indeed, the evolution of the system can be 
described as 
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where G is the conductance of the channel, C 
the compressibility of the reservoirs and \i is the chemical 
potential. The compressibility is analogous to the capac- 
ity of a capacitor. We neglect possible thermoelectric 
effects, since we do not observe a noticeable temperature 
evolution in the reservoirs. 

Because the decay is the slowest process, the deriva- 
tive of the curve around any point is a measurement of 
the current at a certain number difference, where the 
atoms in each reservoir have a known, thermal distribu- 
tion. Therefore, the magnitude of the current measures 
the DC characteristic of the channel. Figure [3j3 shows 
the observed current as function of the number difference 
for the same data set (circles), and for a channel with re- 
duced confinement of 3.2kHz at the center (triangles). A 
linear relation is manifest for both cases, which confirms 
dissipative, ohmic conduction and allows to extract the 
slopes G/C = 2.9(4) s' 1 and 3.7(2) s' 1 respectively. 

The observation of resistance in the ballistic conduc- 
tion shows that the boundary conditions are essential in 
the investigation of transport, as in the Landauer ap- 
proach. Indeed, the free expansion of a non-interacting 



cloud is also ballistic, but the absence of connection to 
reservoirs leads to the absence of any resistance to the 
flow, other than inertia. Furthermore this ballistic ex- 
pansion generally does not depend on the conduction 
properties of the initial cloud, as ballistic expansion has 
even been observed for a band insulator [T2] . 

The Landauer-Biittiker formula states that at zero 
temperature, the conductance of a ballistic conductor 
is equal to 1/h per quantum state contributing to the 
conduction, where h is Planck's constant [20 j. Due to 
the quasi-two-dimensional character of the channel, the 
current is carried by many transverse modes, which are 
not individually resolved because of finite temperature. 
Instead, as the channel confinement is varied, the con- 
ductance is expected to vary linearly with the oscillation 
frequency along the confined direction z, due to the vari- 
ations in the number of modes available. In both mea- 
surements, the reservoirs have the same capacities, thus 
the ratio 0.76(11) of the two slopes is equal to the ratio 
of conductances alone, and agrees qualitatively with the 
ratio 0.82 of trap frequencies along z. We found that the 
linear relation between resistance and trap frequency per- 
sists for various confinement within the accessible range. 
The contact resistance, which naturally appears in the 
Landauer picture, explains the observation of Ohmic con- 
duction even in our defect free channel. While every atom 
that enters the channel on one side exits on the other 
with the same momentum with probability one, only a 
tiny fraction of the atoms from each reservoir can pass 
through the channel at any given time due to the Pauli 
principle. The atoms incident on the channel therefore 
have to redistribute, leading to a drop of chemical poten- 
tial in the contact region, while the chemical potential 
stays constant over the channel [24] [25] . 

To gain further insight into this mechanism, we use 
high-resolution microscopy to observe the density distri- 
bution of atoms in the channel. We do so using in- situ ab- 
sorption imaging along the z-direction, with and without 
current flowing through the channel. A typical picture 
of the density distribution in the channel in the absence 
of current is presented on figure |4]A.. At the sides of the 
picture we observe the contacts with the two reservoirs 
which extend beyond the field of view. Closer to the 
center, the lower column density reveals the presence of 
the channel, which is smoothly connected to the reser- 
voirs. Figure [4^3 shows the difference between two such 
pictures, taken with and without current flowing through 
the channel. We see the small density difference between 
the two reservoirs, which reflects the macroscopic number 
difference shown in figure [2^. 

The red points in figure |4p show the line-density dif- 
ference hi along the channel, obtained by accumulat- 
ing image [4^3 along the x-direction. At the center of 
the channel, the difference is close to zero over 30 /mi, 
while the density difference changes quickly at the sides 
of the channel. This qualitative difference between the 



channel and the contacts indicates a localized chemical 
potential drop, and thus the presence of contact resis- 
tance [26] - Indeed, following the approach of mesoscopic 
physics [25], we introduce a local chemical potential /i(y), 
by requiring that it yields the observed density when used 
in the Fermi-Dirac distribution. As the mean free path 
for atomic collisions is very large, the energy distribu- 
tion of atoms may not be thermal outside the reservoirs. 
Therefore, the definition of \i does not correspond to the 
local density approximation, and does not suppose local 
equilibrium. The variation of the local chemical potential 
is the counterpart of the electric field created by charge 
density variations along an out-of-equilibrium, ballistic 
electric conductor [27] . 

From the current / measured across the channel, we 
deduce from figure pfc the local resistivity [28] 



1 dhi 
nl dy ' 
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where we have introduced the line-compressibility in 
equilibrium K\ = dni/dn, with n\ being the line-density 
along the y-axis. Figure HE presents (in red) the resistiv- 
ity obtained by applying equation (J2| to the in-situ pic- 
ture. The line-compressibility is obtained directly from 
the column density at equilibrium and the shape of the 
trap [28, 29 . At the center of the channel, the local 
resistivity is zero within error bars, whereas it quickly 
increases at the contacts between the channel and the 
reservoirs. We attribute the asymmetry between left and 
right to deviations from an ideal gaussian shape of the 
trapping beam away from the center. 

As the mean free path for atomic collisions is much 
larger than the channel, the resistive current observed 
does not lead to local heating. The energy associated 
with the current is released deep inside the reservoirs. 
The local resistivity observed rather indicates variations 
of momentum, which happen at the length scale of the 
variations of the external potential [25] , leading to a local 
increase of entropy due to the loss of the information on 
the direction of the momentum. 

Many quantities of interest can be extracted from 
the microscopic density distribution. For instance, the 
drift velocity v^ = I/nuj at the center of the channel 
is found to be 200 /ims -1 , or 4 • 10 _3 vp, where vp is 
the Fermi velocity in the reservoirs, which confirms that 
our system realizes the Laudauer paradigm of conduc- 
tion. We also introduce an atomic mobility for the atoms 

VdK ( q d ) , which relates the drift velocity to the 
chemical potential gradient and thus characterizes the 
intrinsic conduction properties of the channel, regardless 
of the density. Figure |4]P presents the atomic mobility as 
obtained from the in-situ pictures for the ballistic chan- 
nel (in red). We observe more than one order of magni- 
tude increase of the mobility at the center of the channel, 
where obtained mobilities diverge, consistent with the ex- 
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FIG. 4: Investigation of conduction using high-resolution microscopy. A: Absorption picture of the density in the channel, for 
a cloud at equilibrium (no current). B: Difference between two pictures taken at equilibrium and with a current of 10 4 s _1 . 
The color is orange for positive difference and green for negative. C: Line-density difference obtained by accumulating B along 
the x-axis, for a ballistic channel (red full circles) and for a diffusive channel having the same conductance (blue open circles). 
The solid lines are smoothed data to guide the eye. D: Focus on the central part of the line density difference. E: Resistivity 
per unit length and time in the channel, as a function of position, for the ballistic (red full circles) and diffusive channel (blue 
open circles), computed from equation (pb. Only the center of the channel is shown, where the extraction procedure for the 
line-compressibility is valid (see Materials and Methods for the details). The solid line indicates zero resistivity. F : Mobility 
in the channel as a function of position for the ballistic (red full circles) and diffusive (blue open circles) channels. The red 
arrows indicate points where the mobility is measured to be infinite. The error bars and shaded region reflect the uncorrelated 
combination of estimated systematic and statistical uncertainties. 



pected infinite mobility of a purely ballistic channel. 

As opposed to the ballistic channel, we have engineered 
a channel where the conduction is diffusive, which is the 
case encountered in typical solid materials. To do so, 
we project a blue-detuned laser speckle pattern onto the 
channel, realizing a quasi-2D disorder [30]. This pat- 
tern has a gaussian envelope with an rms diameter of 
32 /im, an average amplitude of 0.6 jiK at the center and 
a correlation radius of 0.37 /im [28]. We then reduce the 
confinement of the channel down to 1.6 kHz along z, so 
that the atomic conductance of the disordered channel is 
the same as that of the ballistic one studied before. We 
thus have a second system displaying the same macro- 
scopic transport properties, but with a different conduc- 
tion mechanism. The measured line-density difference 
in the disordered channel is shown in blue in figure [4p. 
We observe a continuous decrease from right to left, in- 
dicating a chemical potential drop spread over the whole 
channel, even at the center of the channel as shown in 
[4p. The calculated resistivity, shown in blue in figure^ 
shows moderate variations along the channel, in contrast 
to the ballistic case, and remains sizable at the center. 
The mobility of the disordered channel (in blue in figure 
[4^) is much smaller in the channel than for the ballistic 
case, and clearly saturates to a finite value at the center. 

Our configuration is closely analogous to that of a 



field-effect transistor. The strength of the confinement 
in the channel has been used to vary the conductance by 
changing the density. Further tuning could be obtained 
by adding a repulsive gate laser, tuning the channel to 
fully isolating while keeping an infinite mobility. In ad- 
dition, the effects of disorder in such a device can be 
studied systematically by varying the laser-induced ran- 
dom potential. Metal- insulator transitions, such as two- 
dimensional Anderson localization [31], can be studied 
in a way which is directly analogous to real solid-state 
devices [32]. The ability to further control the disorder 
could be used to study universal conductance fluctua- 
tions |2]. Apart from disorder, various potentials can be 
designed and projected onto the channel using the mi- 
croscope setup [21]. This will allow us to measure the 
conduction properties of various model systems. For ex- 
ample, quantized conduction can be investigated if a sin- 
gle mode can be resolved in the channel [33-35]. Fur- 
thermore, conductance is very sensitive to interactions 
between atoms, and would be an ideal observable to in- 
vestigate strongly correlated fermions. The combination 
of mesoscopic atomic devices with controlled interactions 
opens fascinating perspectives and could shine new light 
on open questions in the field of mesoscopic physics [36J . 
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Materials and Methods 
Potential drop 

Consider a small region of the channel around position 
y. The line density around this point at equilibrium is 
ni(y) and the line density in the presence of current is 
ni(y)+rii(y). 

Since there is no interparticle scattering at the length 
scale of that small region, the energy distribution in the 
presence of a current is not given by the equilibrium 
(Fermi-Dirac) distribution. Indeed, the density varies at 
the short scale of the scattering with impurities or with 
the external potential, while the distribution relaxes to- 
wards equilibrium on the scale of the interparticle inter- 
action mean free path. 

Following the standard procedure in mesoscopic 
physics (see [25] , chap. 2 part 3), we attribute to the re- 
gion around y a chemical potential fi(y), which gives the 
expected density when inserted in the Fermi distribution, 
and coincides with the chemical potential in the reser- 
voirs. In this way, the chemical potential variations are 
local (like the density variations), even if the energy re- 
laxation takes place on a larger length scale, deep in the 
reservoirs. 

We now consider the relation between /i, fii and the 
compressibility k,i = dni/dji. Let f(/n,T,E) be the 
Fermi-Dirac distribution at temperature T, chemical po- 
tential /i and energy E, and g(E) be the density of states, 



m 



= J /(/x, T, E) - /(/i , T, E)g(E)dE (3) 
= (v-»o)^Jf(»o,T,E)g(E)dE (4) 



- Mo) 



dm 



= ^u 



(5) 
(6) 



where we have introduced the local chemical potential 
at equilibrium /i > which follows from the local density 
approximation at equilibrium, and jl = /i — /io the local 
deviation of the chemical potential away from equilib- 
rium. 

Knowing the current / flowing through the channel, we 



p{y) 



(7) 



now introduce a local resistivity 

djl 1 
~dyT 

which reflects the local scattering with obstacles. Al- 
though the Joule heating happens deep in the reservoirs 
where excited atoms release their energy via collisions 
with other atoms, the momentum randomization due to 
scattering with impurities is local, and leads to a local 
increase of entropy as information on the direction of 
motion of particles gets lost. 

Experimentally, to get the derivative dhi/dy we fit a 
line to the line density difference at each position y within 
a window of ±9 /mi around this point. 

To extract k>i at position y we make use of the isotropy 
of the compressibility. We consider at fixed position y 
the variation of n co i along the x direction and compare it 
with the known trap shape [29 . Consider the gas around 
position y: 



Ki(y) = 



dni(y) 



= g- n(n(x,y,z))dxdz 
= — {x,y,z)dxdz 

J J dx \dx J 

dn fdV\~ x , , 
dxdz 



dx V dx 



(8) 

(9) 

(10) 

(11) 
(12) 



The confinement along the x direction is ensured by the 
optical dipole trap, and is constant over the channel. 
Close to the center of the channel, the confinement along 
z is mainly ensured by the two repulsive lobes of the laser 
beam. Therefore, the variations of the potential along the 
z direction is independent of the variations along the x 
direction and can therefore be pulled out of the z inte- 
gration. 



Ki(y) 



dV 
dx 

dVy 1 dn co i 
dx J dx 



T^dz ) dx 
dx 



dx , 



(13) 
(14) 



where n co i is the column density. 

A slice of ±9 /mi around one position y is taken from 
the density picture at equilibrium (figure 4A) and the 
average along the y direction is calculated. This results 
in a one-dimensional density profile along the x direction. 
We fit a gaussian to such a profile and use it in equation 



(14), in order to avoid the noise generated by numerical 
differentiations and ratios of those 1291. This fits the 



shape of the cloud within the errorbars. In equation ( 14 ) , 



V is taken to be the known gaussian shaped trapping 
potential. 



Mobility 

In electric conduction the mobility relates the drift ve- 
locity Vd to a potential gradient (electric field). In our 
case the gradient in the chemical potential plays the role 
of the potential gradient and the atomic mobility is 



Vd- 



dy 



Vd • Kl 



dni 
dy 



(15) 



where we used equation (pi). 

In the data analysis we attribute an infinite atomic 
mobility to points where ^ is zero or slightly negative. 
These points are indicated with red arrows in figure 4. 
The upper limit of the red shaded region is as well set to 
infinity if the lower edge of the error on -^r extents to 



negative values. 



dy 



Properties of the confining potential 

The channel is imprinted on the atoms using a 
532 nm wavelength laser beam, with waists 30.2(3) and 
10.3(3) /im along the y and z directions, respectively. 
This beam passes through a holographic plate (Silios) 
dephasing the upper part of the beam by tt with respect 
to the lower part [22, 23 . The beam has been character- 
ized in a test setup and shows a slight asymmetry along 
the y axis away from the center. The contrast of the cen- 
tral region of the beam compared to the lobes is larger 
than 0.99. In the parameter regime explored in the pa- 
per, the residual light leads to a repulsive potential along 
the y-axis smaller than 40 nK, much smaller than the os- 
cillation frequency along z. 

The oscillation frequency along the tightly confining 
direction z has been measured in-situ using parametric 
heating on a microscopic cloud [21], and the measured 
frequency agrees with the fitted curvature of the intensity 
profile observed on the test setup. Along the propaga- 
tion direction of the beam, the measured curvature of the 
intensity profile varies by less than 5% on a length scale 
of 200 //m, and is therefore constant over the transverse 
radius of the cloud. 



Properties of the disordered potential 

The disordered potential is created by passing through 
a light diffuser (Luminit) with a beam at 532 nm, then 
through a high resolution microscope objective [21], cor- 
rected for aberrations at that wavelength. The resulting 
disordered potential is observed directly in-situ using a 
second microscope identical to the first. In this way, we 
directly characterize the potential correlation properties 
by taking pictures of the potential as seen by the atoms in 



the glass cell. Figure SI presents a typical observation of 
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FIG. SI: Properties of the disordered potential. A : central region of the disordered potential as observed with the high- 
resolution microscope used to image the atoms. The total width of the picture is 30 /xm. B : correlation function of the 
intensity distribution observed in A. C : horizontal cut of B, showing the correlation properties of the disordered potential (in 
blue). The solid green line shows a Gaussian fit yielding a measured correlation radius of 460 nm. 



the potential characteristics. Figure |Sl} \ shows a zoom 
on the central part of the disordered potential where cor- 
relation properties are computed. Figure [Slp> shows the 
measured correlation function obtained using the inverse 
Fourier transform of the power spectrum of A. Figure 



Sip shows a cut along the horizontal axis, of figure B, 
together with a Gaussian fit, yielding an observed corre- 
lation radius (1/y/e radius) in this direction of 460 nm. In 
the other direction, the same fit yields an observed corre- 
lation radius of 560 nm. To obtain a faithful estimate of 
the correlation properties of the disorder on the atoms, 
we deconvolve the correlation properties using the mea- 
sured, gaussian point spread function of the microscope 
used to image the pattern [2T] , This yields a correlation 
radius averaged for the two directions of 370 nm. In ad- 
dition, using the in-situ images of the light intensity, we 
also observed the expected exponential probability dis- 
tribution of intensities. 

The envelope of the disorder potential is straightfor- 
wardly obtained from a gaussian fit of the observed pro- 
file. The total power sent through the microscope, to- 
gether with the average fitted profile of the envelope, 
yields the average depth of the disordered potential cited 
in the text. 



